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SUL'IMARY- 



The differential equation of Chaplygin^s jet problem 
is utilized to give a systematic development of particular 
soluticx^s of the hodograph flow equations, which extends 
the treatment of Chaplygln into the supersonic range and 
completes the set of particular solutions. 

The particular solutions serve to place on a rea- 
sonable basis the use of velocity correction formulas for 
the com.parison of incompressible and compressible flows. 
It is shown that the geometric-mean type of velocity 
correction formula introduced in an earlier paper, part I, 
has significance as an over-all type of approximation in 
the subsonic range, 

A brief review of general conditions limiting the 
potential flow of an adiabatic compressible fluid is 
given end application is made to the particular solutions, 
yielding conditions for the existence of singular loci 
in the supersonic range. 

The combining of particular solutions in accordance 
with prescribed boundary- flow conditions is not treated 
in the present paper. 



INTRODUCTION 



This paper presents a theoretical investigation that 
may be regarded as a continuation of studies initiated 
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In part I (reference 1). In part I an attempt was made 
to unify the results of Ghaplygin, von Karman and Tsien, 
Temple and Ysrwood, and Prandtl and Glauert insofar as 
their results were concerned with velocity and pressure 
correction factors for the correspondence of incompres- 
sible and compressible flows. In addition, two new 
velocity correction formulas were introduced that appeared 
to have a somewhat wider range of applicability than the 
formulas of the afore-mentioned authors. Most of the 
results of part I were obtained with the use of two par- 
ticular solutions of the hodograph equations. These tv/o 
basic solutions correspond to a vortex and a source in a 
ccmpre s s ible fluid • 

It was mentioned in part r that, in order to treat 

the exact boundary problem of uniform flow of a compres- 
sible fluid past a prescribed body, a general set of 
particular solutions of the hodograph equations had to be 
obtained. Such a study is given in the present paper, 
which incidentally helps to clarify the nature of the 
velocity correction factors of part I - in particular, 
the one referred to as the ^^geom.etric-mean'^ type of 
approximation. In addition, many interesting types of 
flows are disclosed from a physical interpretation of 
the particular solutions. A few such solutions hsve 
already been obtained and discussed by Ringleb (refer- 
ence 2 ) . 

Several mathematical approaches exist by means of 
which -particular integrals of the hodograph equations 
may be obtained. Two such approaches, mentioned in 
part I, m^ay be attributed to Chaplygin (reference 5) and 
Bers and Gelbart (reference I4.) and are analogous to an 
exponential and to a power-series approach, respectively. 
Another method of defining particular integrals is the 
Integral-operator method of Bergman, (reference 5)- In 
the present paper the differential equation,- first used 
by Chaplygin in his treatment of jets (reference 3)* 
provides the bssis for the definition of a complete set 
of particular solutions. 

The scope of the present paper lis limited chiefly 
to a systematic study of the fundamental solutions and 
to the physical interpretation of some of the particular 
flows represented by them. The combining of particular 
solutions to represent uniform flow past a prescribed 
body is not treated herein. It is believed, however, 
that the present study may serve as a basis for further 
development and clarification of this important problem. 
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s™bols 

X, y rectangular coordinetbs in plane of flow 

q magnitude of fluid velocity 

e angle included by velocity vector and posi- 

tive direction of x-axis 

p density of fluid 

p pressure in fluid 

a velocity of sound in fluid 

M Mach number (q/a) 

Pc^ Po> ^o quantities referred to stagnation point q = 0 

^ velocity potential 

V stream function 

Y ratio of specific heats (approx. l.ii for air) 
1 

P - — - — (approx. 5/2 for air) 
Y - 1 

2paQ maximum fluid velocity (corresponding to 
P = p = a = 0) 

dlmensionless speed variable 



Tg sonic value of T (t = ; approx. 

1/6 for airj 

?or p>0 (or 1<Y<*)> the range of t is 

0 ^ T ^ 1. ■ 



M 
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GBMSRAL - PARTICULAR SOLUTIONS OF TEE HODOGRAPH EQUATIONS 

Hodcgraph Equations 



The linear equations in the hodogrsph V8ri8b].es 8 
and q, which relate the velocity potential 0 and the 
stream function \j/ for the steady two-dimensional flov/ 
of a nonviscous compressible fliiid, are 



> 



(1) 



in which, for the adiebr.tic equation of state between 
pressure and density, 



^ q 



(1 - T)' 



and 



ciqvpqy 
- Po /-I _ f,2\ 

""q(l.T)P^^ 

(See equations (21) and (25) of reference 1,) 

In the Incompressible case t — > 0^ equations (1) 
can be expressed in the Cauchy-Riemann form. Pe.rticule.r 
solutions Q - 0 + iw can be' expressed in this case es 
any analytic function of the complex vsriable 

. ■ w = 8 + i log q (2) 
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or as any analytic function of the related exponential 

function 



Thus, an infinite set of particular Integrals of equa- 
tions (1), in the incompressible case, referred to herein 
as "the powers set,'' is w^^. ^/^/hen k is a positive 
integer, the particular solixtions vanish at the origin 
(9 =0, log q = 0) and, when k is a negative integer, 
the particular solutions are infinite at the origin. In 
the case of nonintegral values of k, the origin is a 
branch point of the functions w^'^. 

Another infinite set of particular integrals of 
equations (1) in the incorapressible case, referred to 
herein as "the exponential set,'- is 



where, again, k can take on any value - integral, non- 
integral, positive, or negative. 

In the compressible case, the particular solutions 
corresponding to the powers set w^^^ (that is, the par- 
ticular solutions which reduce to w^^^ in the incompres- 
sible case T — >0) depend on v/hether the coefficient of 

is real or imaginary - a consequence of the fact 
that, in the compressible' c^se, 0 and >J/ do not 
satisfy the same differential equationc For example, 
for k 1, the two functions corresponding to w 
and iw, which have been developed in part I, are 




W = e + iL 



and 




v/here 



T 



log q + f(T) 



and 



L = log q + g(T) 



6 



NACA ARR No, 11+129 



and f(T) and. g(T) each vanish for t = 0. (See 
equations {26) and (27) of reference 1.) 

The development of other functions corresponding to 
the pov^er set w^^, for positive integral values of k, 
follows according to the method of Bars and Gelbart. 
(See expression (22) of reference 1.) Since the present 
paper is chiefly concerned with the ^ functions corre- 
sponding to the exponential set e"^^^'^ , the powers set 
is not further discussed. 



Chaplygln Differential Equation 

The Functions P^^ and Qj^ 

Corresponding to the exponential sets in the incom.- 
pressible case 

e"^^^ = cos k8 - iq^ sin k9 



and 



le"^^^* = q^ sin k8 + iq''^ cos kG 



there appear in the compressible case functions designated, 
respectively, 

Pl^(q) cos ke - 1 Q^(q) sin kG 

and 

Pj^(q) sin ke + 1 ^i^iq) cos k8 

where the functions Pi^(q) and ^^^(q) satisfy second- 
order differential equations. These equations are easily 
obtained by substituting in equations (1) the product- 
type solutions 



4 = ?k(^) sin (1^9) 



= Q, (q) It^' (.-ke) 

k k cos 



> (1+) 
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In view of equations (1) it is observed that 



k P,(q) = Eo , 



dPk(q) 



dq 



dq 



> 



(5) 



The functions Qj^(q) satisfy the second-order differ- 
ential equation 



dq^ 



+ (l + .M^)q ^ - k^(l - M^) = 0 (6) 



The functions P]^(q) can be obtained from Qij(q) by 
means of the first of equations (5). Equation (6) may 
be reduced to a standard type by introducing t as the 
independent variable. Put 



= q-" Yi,(t) 



(7) 



where clearly Y^{t) — >1 as t — >0 (incompressible 
case). With the use of the symbolic relations 



q d = 2t j_ 
dq dT 



„2 _ d-^ 
dq^ dT^ 



+ 2t 



dT 



and the relation 



M' 



2 _ 2p- 



1 - T 



the desired differential equation is 
p [(k+1) -(k+1- ,«)Tj 

dT' 



T(l-T) [(k+l)-(k+l-,«)T] -^+|-pk(k+l)\»0 



(8) 



8 
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Equation (S), which is of the hype rgeome trie tjpe, was 
first introduced by Cheplygin in his memoir on gas jetii 
(reference 5), 



The Functions Yj^ and Y.j^ 

Chaplygin treated the subsonic flow of a compres- 
sible fluid through jets with straight-line boundaries. 
For such problems the hodograph variables 9 and q 
are natural variables in the sense that the solid and 
fluid boundaries are described by 6 = Constant and 
q = Constant, respectively, and only the particular solu' 
tions of equation (8) with positive characteristic 
index k are needed. In the present paper a complete 
ordered set of particular solutions of equation (8) is 
obtained, which extends the results of Chaplygin into 
the supersonic range and ^ to negative values of the 
index k. Two types of solutions of equation (8) for 
nonintegral values of k are 

Y^rir) = F(ak, b^^, k+1; r) (9) 



ana 



where 

+ bj, = k - 



and 



?(a, b, cj t) = 1 + i2.T+ ' • + 



a(a + l)b(b + 1) P 



0 



21 c(c + 1) 



It is now shown that only one of the solutions need be 
used. For positive values of k, the requirement that 
Y',^(0) = 1 excludes the use of equation JlO ) . For nega- 
tive values of the index, the solution Q'^i^q) = q" Y^j^{ t) 
obtained with the aid of equation (10) is,'" except for a 
constant factor, equivalent to the solution ^(q.}--^q^ ^k^'^'^ 
obtained with the aid of equation (9). Thus 
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= q"^T-^ F(a.i^+k, b.^^+k, k+1; t) 
= q-^T^^ p(a^., b-^, k+1; t) 



Ilenoe, only the solutions given by equation (9) ai'e 
needed for the detennlnetion of Qj^(q) end Q.]^(q)- 
Then 

= q^ F(a^, bj^, k+1: t) (11) 

and 

= q"^ ^Kk' ^-k^ --^^^5 0 

= q-^ P(ak-k, b^-k, -k+l; t) (12) 

Observe that both types of hypergeometrlc functions 
appearing in equations (9) and (10) are utilized in the 
expressions for Qi^(q) snd Q^]^(q). 

The foregoing discussion has been limited to non- 
integral values of the index, positive or negative. When 
the index is integral and positive, equations (9) and (11) 
remain valid. When the index is Integral snd negative, 
however, equation (12) does not in general lead to a 
meaningful solution and consequently another independent 
solution is to be sought. The desired solution for 
Y_j^(t) in such cases contains a logarithmic term and 
again is subject to the condition that it reduce to unity 
for T = 0 (incompressible case). The expression for 
Q_^(q) is then given by 

Q.k(q) = q-^' Y.k(^) ' (^5) 
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Cese of Y " 

Consider as an example the von Karman-Tsien treat 
ment of compressible flow (reference 6) in which the 

adiabatic index y = -1 or p = -^^^ Then 



k + 1 



b,, = ^ 



For a negative integral index, equation (15) may appear 

to be applicable, in which case the expression for 
Y_-^(t) would be a polynomial of degree k - 1. An 

examination of equation (12) shows ^ however, that for this 
case no infinities arise and that, when the index is 
negative, integral or nonintegral. 



The hypergeometric series represented by Yi.(t) con- 
verges for values 0 = jil < 1. For the present case of 
Y ~ -1 or p = values of t corresponding to 

positive values of M lie outside the range of conver- 
gence, A closed expression for Y ^^{f) can be found, 
however, . for this case which, by analytic continuation, 
is therefore valid for all values of " t. Thus 

Y_i^(T) = F(i--i^, -|, l-kj t) 

k 



Similarly, froin equation (9), when the index is positive. 



1 + (1- T)^^ 
2 
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Observe that 



2 ' 2 
1 + (1 - T)V2 



-k 



(q) 



-ik 



1 ^ (1 - T)^/l 



This Identity for the von Karman-Tslen case corresponds 
to the identity = -i— for the incompressible case. 

Case of k = 1: 
For k = 1, 

a-|^ = 1 = -p c-]_ = 2 

Then, for the positive index, 



= q P(l, -p, 2; T) 

1 - (1 ^ T)P-^^ 
= q i 

(P + 1)t 



(15) 



For the negative integral index, it may appear st first 
glance that equation (1$) is needed; however, equation (12) 
does yield a relevant and finite result and accordingly 
is the equation to be used. Thus 
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Lim F^a^r-k, b^-k, 1-k; t) = 1 + i t ^ — + ^^^^ " "^^ 

k->l ^ ^ 2 2 X 2J 2x3.' 



2 X hi 

= i.i-^ri-(i-T)p^^l 

2 S + 1 L J 



and therefore 



Case of k = 0: 

The exceptional case of k = 0 is directly treated 
by means of equation (8). The differential equation for 
Yq(t ) or Qo(t ) then is 



d_ 
dT 



(1 - T)P dT 



= 0 



The general solution of this equation can be written as 



Qo(q) = 2C^ log q + Ci / [(1 - t)^ - i] ^ + 

Jo 
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'.vhere C-^ and C2 are arbitrary constents of integra- 
tion. The constents C]_ and C2 are determined by the 
imposed condition that the expression for Qq^*^^ reduce 
in the incompressible case simply to log q. Then 

Cl - 2 



end therefore 



C2 = 0 



Qo(q) = log q + I r[(l - T)P - 1] ^ 



(17) 



In a similar manner, from the differential equation for 
Po>- 



_d_ 
dT 



T(l - T)'^"'^ dPo 

1 - (2p + 1)t dT 



= 0 



the expression for P^- is obtained as 



Pp(q) = log q 



vJo 



1 (2(3 + 1)t _ ^ 



a - T) 



p+1 



dT 



(18) 



It is remarked that the functions Qo(q) and Po(q) 
are identical with the elementary functions L(q) 
and L(q), respectively, introduced in part I (refer- 
ence 1) and are associated v/ith a vortex and a source 
type of flow. 



The Functions Rj^ and Sj^ 

A linear homogeneous differential equation of 
order n can, in general, be reduced to a differential 
equation of order n - 1 by means of an exponential- 
type substitution for the dependent variable. Chaplygin 
made use of such a substitution to reduce the second- 
order differential equations satisfied by P, and 
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to first-order equations of the Riccati fom, In order 
to study properties of the functions P-^ and in the 

subsonic range for only positive values of k. In the 
present analysis the Riccati equations are also found 
useful in order to extend the study of the functions Pj^ 
and to the supersonic range for both positive and 

negative values of the index k. 



The second-order differential equations for 
and Q^, with t as the independent variable, are 



and 



_d_ 
dT 



t(1 - T)P+1 dP 



1 - (2(3 + 1)t 



dT 



k2 (1 ^ t)P 
h 



= 0 



_d_ 

dT 



T 

(1 - T ) P • dT 



. kl 1 - (2P -t- l)T 



T(l - T) 



The corresponding first-order Riccati equations are 
obtained by substituting for P,^ and new dependent 

variables Yi-^ and S^., respectively, es follov/s: 



RkdT 



or 



Rk = -r- 



2t 1 dP, 



k Pjj. dT 



2t d T 

— — log Pj^ 

k dT ^ 



(19) 



and 



Qk 
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or 



^ k Qj^ dT 
- 2t d T ^ 



(20) 



The equations satisfied by Rj^Ct) and Sj^(t) are 



dP.i 



^ 1 + (2p + l)T p 



dT 1 - (2p + 1)T 1 - T 
and 



Ri 



2 1-(2(3^-1)t 
^ " 1 - T 



= 0 



(21) 



dSk 
dT 



1 - T 



2t 



2 1 - (2|3 + 1)t 



i - t 



= 0 



(22) 



Initial conditions for Ri^(T ) and Sj^(t) are found by 
examinstion of the incompressible case t — >0. In this 

case Pk = Qk - since 2t ^ = q it follows 

from equations (I9) and (20) that 



R^(0) = Si,(0)-= 1 



The following important relation exists between the 
functions Rk(T) and Sj^(t): 

^ ^ 1 - T 



= 1 - (25) 



Equation (2$) can be verified directly from the hodograph 
equations (1). It may be noted at this point that this 
result Is of significance in connection with the 
geometric -me an type of velocity correction factor intro- 
duced in pert I and is discussed more fully in a later 
section. 
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Before the functions R]^{r) and S^{r) are treated, 
certain general observations can be made regarding the 
functions V^^ir), Q^^ir), end S^{t). Chaplygin, 

who limited his investigations to the subsonic range and 
to positive values of the index has shown that 

and consequently the other functions possess no roots 
for any value of the independent variable in the subsonic 
range, with M - 0 excluded. In the supersonic range 
M > 1, P^^t) and Q],^(t) in general possess zeros. 
Certain relations obtained by means of equations (19), 
(20), and (23) between p,^, q^, r^, and at the 

zeros of Pj^ and . are sujnmsrized as follows: 







2lk 
dT 


dT 






0 


Max or min 


• • • 


0 


CO 


0 


Max or min. 


0 


0 


• c • 


0 


CO 



It is remarked thst the number of zeros of Q^, 
83 a function of the index k, can be found from an 
expression developed by Klein and Hurwitz (reference 7) 
in connection with the zeros of the hype rgeome trie func- 
tion. In general, the number of zeros increases with 
the magnitude of the index k and is inl'^inite for 
k = ±00, 

A further observation of interest can be made in con- 
nection with equation (2^). Chaplygin has shown that, 
for positive finite values of k (and the same is true 
for negative finite values of k), the functions S^(t) 
are not zero for the sonic value t = Tg or M = 1 . 
Prom equation (23) then, it follows that the functions 
R^(t) = 0 for M = 1. 

In view of the relation between the functions R^^ 
and given by equation (23), only Sj^ need be dis- 

cussed. The Riccati equation (22) may be" used to discuss 
certain properties of the fujiction S]^ but in general, 
for numerical evalurtion, the original definition (equa- 
tion (20)) in terms of the function q^^ may be used 
directly; ^ 
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2r_ 

k Qj^ dT 



or 



Yk dT 



k 



In general, the functions S]^ are expressible in InTl- 
nite series. For* several values of k, however, S 
can be expressed in closed forms. For k = 0 and 
k = t^, S]^ may be obtained by a limiting process from 
equation (20 )j however, for these special cases the 
Riccati equation (equation (22)) yields the results 
directly. Thus 



So = (1 - T)P 



(21+) 



and 



>±00 



1 ~ (2g + 1)t 



1 - T 



1/2 



(25) 



The cases k = 1 and k = -1 may also be expressed in 
closed .form. With the aid of the equations (I5) and (I6) 
for Q-^ and Q^-^, equation (20) yields 



= 1 - 2 



1 - (1 + pT)(l - T) 



1 - (1-T) 



P+1 



(26) 



and 



3., = 1 <2„ 

2 p + 1^ -J 

In order to illustrate the behavior of some of the 
functions thus far introduced, a n^araber of tables and 
figures are given. All the calculations heve been 
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performed with the adi^ibatic index ^ - Table 1 

gives velues of y^. as a function of M or t for 
several positive and negative values of the index k. 

Figure 1 shows the Yj. functions plotted against M. 

Values of the S^^ anci R^^ functions are given in 

tables 2 and 3 ?aid are plotted against M in figures 2 
and $. 



The Functions fj^(T) and gj^(T) 

In the incompressible case, the sets of functions 
and Py^ can be reduced to a single function log q by 

means of a simple operetor ^ log. Thus 

"k ^k - 

and 

1 . k 

log q^^ = lo 



a- 



q 



Thi3 seme operation applied to the functions Q^^ and 
in the compressible esse serves to define two u'seful sets 
of functions log q + t^{r) and log q + gi^(T), respec- 
tively. Thus 

^ log = log q + f^(T) (28) 

and 

~ log = log q + gj^(T) (29) 
Prom equation (7), namely. 



it follows that 



f^{r) = I log Yy{T) (50) 
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Prom equation (5) for Pj^. . and equation (20), which 
defines S^, 



(1 - T)P ^ 



It follows that 

gk(T) = f log 



1 , Y^(t) S^(t) 



f^ir) + - log 



(1 - T)' 



(51) 



For exsmple, for k = 1 and k = -1 and with the use of 
equations (15) and (16), 



r.ir) = log 



1 ~ (1 - T)P 

(P + 1)t 



+1 



(1 - T)P[1 + (2(3 + 1)t] - 1 
S-^ir) = log 

(P + 1)T(1 - T)P. 
f_-L(T) = - l0g<jl + [l - (1 - T)P-'1] 



(52) 



(55) 



(54) 



g_]_( t) = - log 



(5P + 1) - P(l + T)(l - T)P 
2(p + 1)(1 - T)P 



(55) 



For k = 0 and k = ±<», equations (50) and (51) require 
a limiting process for their evaluation. Alternate forms 
for fk(T ) and S\ii'^ ) t)e obtained, hov/ever, by 

means of equations (19) and (20) defining Rif(''') and 
Sj^(t), which yield the results for k = 0 and k = 
directly. Thus . 
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and 



Si^(t) = 



Jo 



1 rr 



2 / L^^k(^) 



- 1 



" dT 



J T 



(36) 



(57) 



where Ri^(T) and Sj^(t) are related according to equa- 
tion (25). Then 



(36) 



1 - <2fi * 1)T 

(1 .t yP-^i 



dT 

T 



(59) 



and 



f±oj(T) = £iod(T) =^ J 



1 1 



1 - T 



1/2 



T 



(I|.0) 



It is v/orthy of special notice that the functions ^o{t), 
g^ir), and %^(t) are identical with the functions 
-(t), g(T), and h{T), respectively, which formed the 
basis of part I (reference 1). In addition, the expres- 
sions log q + fo(T), log q + SqCt), ^and log ^^f^i^) 
are Identical with the functions L, L, and H, 
respectively, which were introduced in part I. 

A number of functions and g^^ have been 

calculated, with y ^ ^or several positive and 

negative values of the index k, and the values are 
given in tables i4. and 5 snd plotted in figures ij. and 5. 
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The opportunity is taken here to note that, for the 
von Ka'rma'n-Tsien case 



(y = -1 or p = -|), 

.1/2 



i'k = Sk = - los 



1 + (1 - T)- 



= log. 



1 + 



and that the sets of functions and Qi^, as in the 

incompressible case, are reduced to a single function by 

the operator -log; namely (compare equation (11+)), 
k 

2yin? 



log q + log 



1 + vi - ^ 



M 



,2 



In fact, the complex -flow' potential ^ + it can be 
expressed as an anelytic f unct ion of a single complex 

variable 9 + i log q V - ^- Tsien has made use 

1 + vi - 

of this complex variable in his hodograph treatment of 
the compressible flow past an elliptic cylinder (refer- 
ence 8) . 

Velocity Correction Factor 

The solution of the problem of an exact corre- 
snondence between the flow past a prescribed body in an 
incompressible fluid and the flow psst the same body in 
a compressible fluid is a difficult matter. This 
problem can be solved exactly for certain types of flow 
patterns- (not past closed shapes), such as flows inside 
or outside angles or channels, and for certain flow 
. singularities such as a vortex, source, and doublet - 
types of flow which can be associated with the particular 
solutions G^. Some of these types of flow are illus- 
trated by examples in the following section. Combining 
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particular solutions to represent uniform flow past a 
prescribed body is a complicated process, since^ the treat- 
ment of infinite series in the fvmctions for both 
positive and negative values of k is involved. Further- 
more, the process of - returning to the physical-plane 
variables from the hodograph-plane variables hinges on 
nonelementary parts of differential geom.etry. Certain 
types of jet problems can be properly treated in the sub- 
sonic range by series in with k positive, as was 
shown by Chaplygin (reference 3). Thus, it appears that 
much work remains to be done in order to render feasible 
exact and practical solutions for uniform flow past pre- 
scribed bodies in a compressible fluid. Because of the 
difficulty and complexity of the general problem of flow 
in a compressible fluid, attempts have been made by a 
number of investigators to obtain results by means of • 
velocity correction formulas that serve to place in 
correspondence velocities in an incompressible and in a 
compressible fluid. 

In part I the velocity correction factor was dis- 
cussed with particular reference to the two functions L 
and L (Qq and of the present paper) associated 

with a vortex and source type of flow, respectively. The 
main justification for the results of part I was the 
yielding ^and the unifying of the results of Chaplygin, 
von Karman and Tsien, Temple and Yarwood, and Prandtl 
and Glauert. The knowledge of the infinite set of func- 
tions ?ir and discussed in the present paper can 
now serve to establish further on a reasonable basis the 
concept of a velocity correction fomula. 

In order that a single velocity correction factor 
be feasible, even for a flow associated with a particular 
solution, it is necessary that P^^ s Qj^. Consider, for 
example, the functions Qj^ and P^ Insofar as the first 
power of the variable t is concerned. It can be shown 
easily that 



- log Qj^ = log q + f^Cr) 



21+ 
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and 



i. log = log q + gj^(T ) 



- log q - jpr 



Thus, to the first power of t and independent of k. 



fw(T) = gl,(T) 



4pt 



Then 



= Qk 



( -l-T 



The nature of the correspondence between the Incompressible 
flow and the compressible flow is such that 



(ill) 



Without going into any details here of the field point 
correspondence or of the boundary distortion, the veloci- 
ties in the incompressible and compressible cases may be 
placed in correspondence as follows: 



(log q)^ = (los, q - |-PT ) 



or 



Qi = %^ 



"2 Pt 



(il2) 
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This result implies that the complex variable 

6 + 1 (log q - 2'^'^} the compressible case corresponds 

to the complex variable 6 + i log q in the incompres- 
sible case. Equation (I|.2) represents the approximation 
of Temple and Yarwood discussed in part I. 

Consider now the functions Qj^ and insofar as 

large values of the index k are concerned. It is 
recalled that, as the index k — > ±<», 

and that 

1 log 1 lOS 

— > log q + h(T) 

where 



h(T) = f^Jr) =; g^^(T) 



Then, as k — > ±00, 



k 



The function h(T) is expressed in integral form in 
equation ([{.O) and has been evaluated and tabulated in 
part I. (See also table 1^ and figo l\. of the present 
paper.) The correspondence of velocities in the incom- 
pressible and the compressible case is given by 



^i = 



h(T) 



(i|3) 
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Equation (^5) constitutes the geometric -me an velocity 
correction formula introduced in part I and is limited 
to the subsonic range 0 ^ M ^ 1. It is observed that, 
for positive values of k, h(T) lies betv/een 

and Si^i'^) in magnitude. Moreover, the deviation of 

^h(T) fjoQm e'^^^'^^ and e^^^'^^ is quite small in the 
subsonic range. (See table 6.) 

The foregoing remarks, together with the fact th^t 
the geometric-mean type of approximation contains the 
results of Chaplygin, von Karman and Tsien, Temple and 
Yarwood, and in the limiting case of small disturbances 
to the main flow the exact Pr pndtl-G-lauert rule, lead.- 
to the suggestion that it may be adopted as an over-all 
type of approximation in the subsonic range. 



Flow Patterns Corresponding to the Particular Solutions 

Before the flow patterns corresponding to the par- 
ticular solutions '^^k given by equations 
for the compressible case are discussed, it is instruc- 
tive to examine the incompressible case. Consider the 
complex velocity potential 



Q = -f i\|/ 



= Uz 



n 



ihk) 



where U and n are constants and z = x + iy^ It is 
well knovm that, if ^ ~ ^ v/here a is an sngle between 

0 and 2Tr, equation {hh) represents the flow in a sharp 
angle. For example, the flow inside a right angle is 
obtained with n = 2 and the flovv outside a right angle 

p 

is obtained with n = f-. Again, the value n = 1 or 

a = n corresponds to a uniform flow and the .value n = 2 
or a = 2Tr corresponds to the flow around a semi- 
infinite line. Clearly, all the angle flows are obtained 
Y\^ith values of n between l/2 and oo. Other types of 
flows are given by other values of n. For example, 
n = -1 corresponds to a doublet and the remaining 
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iiegatlve integers are assocleted with singularities of 
higher order than the doublet. In addition to the flows 
described by the powers z^, there are the two funda- 
mentel flows, the source and the vortex, associated with 
the function log z. If, now, it is desired to obtain 
generalizations for the compressible case of the fore- 
going particular flows, the procedure is first to express 
52 or a)/ for the incompressible flov\^ as a function of 
the hodograph variables q and 6 and then to replace 
q 'oy Pk or Q^, respectively. Several examples will 
best illustrate this procedure: 

(1) Consider the compr'essible generalization of the 
angle flows. By means of the relation 

dQ -iw 
~ = e 

dz 

where w = 9 + 1 log q, the hodograph complex variable w 
is Introduced as independent variable in pl&ce of z. 
Prom equation (ijij.) 



dz 

= e-i^ 



Hence 



1 

n-1 



(Un) 



n-1 



and 
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Then 



n 



n 




n-1 



(Un) 



If Y is replaced by k, the compressible general! 
z at ion of the angle flov/s is given by 



\ = - 



k 



Qt, sin k6 



(1+5) 




The inside angle flows are given by values of k in the 
range 1 < k < «> and the outside angle flows, by values 
of k in the range 1 ^ -k < For example, k = 2 

for the ilovv inside a right angle, and k = -2 for the 
flow outside a right angle. Other types of flow are 
given by values of k in the range -1 < k ^ 1. 

The case k = 1 or n =. ±«> is exceptional and, in 
fact, corresponds to the incompressible flow 



vi^here c is a constant, 

(2) Consider the compressible generalization of the 
doublet. The complex velocity potential for the incom- 
pressible doublet at the origin is 



Q = e 



cz 



(1+6) 



z 



The ref lected-velocity vector is 



dO 
dz 



- e 
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Kence 

z = le"^ 



and 

^1 

Q = -3.e 



The stre8m function for the incompressible doublet is 
then given by 

1 

\|/ = - cos 3-9 



The compressible generalization of the doublet is 
therefore 

(5) Consider the compressible generalizetion of the 
source. The complex velocity potential for a unit source 
at the origin is 

Q = log z 

The ref lected-velocity vector is 

dz ^ 

- ^-iw 

— e 

Hence 

z = eiw 



and 



Q = Iw 
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The velocity potential for the incompressible source is 

= - log q 

The compressible generalization of the source is then 
given by 




(k) Consider the compressible generalization of a 
point vortex. The complex velocity potential for a 
vortex of unit strength at the origin is 

0 = - i log z 

The reflected- velocity vector is 

dQ ^ ^ i 

Kence 

z = -iei'^ 



and, except for an additive constant, 

n = w 

The stream function for the incompressible vortex is 

^ = log q 

The compressible generalization of the vortex is then 

given by 
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Trsnsf ormation from the Hodograph 

to the Phy s ic al Var 1 abl e s 

Given the velocity potential ^ and the stream 
function in terms of the hodogreph variables 8 

and q, it is possible to express the ccordlnstes x 
and y of the physical plane in terms of 9 and q. 

Prom the basic flow equations 

6x p 6y 



6^ ^ ^ Po 
6y p 6x 

it follows (see equation (6) of reference 1) that 



The real and imaginary parts of this equation yield 

, sin 9 ) do 
P oq 



q 



V6q 



^ cos 9 - ^ 1^ sin e) dq 



+ cos 9 - 'n sin e) de 



V69 



p 66 



> 



(47) 



^ (H « * ^ If 0 



d9 



Equations ([[.7) relate the differential line elements in 
the physical x, y plane and the hodograph 9, q plane. 
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vvhen expressions for ^ and \j/ as functions of 9 
and q are known for a given flow^ the integrals of 
equations (ij.7) are the equations of transformation of 
the 9, q to the x, y coordinates. It may be remarked 
that the hodograph flov^^ equations (1) are the inte^ra- 
bilitjr conditions for the differential equations (1+7). 
The right-hand sides of equations (i|7) are therefore 
perfect differentials. 

Consider one set of particular solutions from equa- 
tions 

^ .= - Qk(q) sin k6 

where k = ±1 and k = 0 are , excluded. By the use of 
equations (5), it can easily be verified that 



k 



Pk Po 0 ) co3(k -»• 1)6 

,q ' pq V k + 1 



cos(k - X)9 
k - 1 



!k . P^q/\ sln(k + 1)9 




+ Constant 



pq 



k 



k + 1 



P 

pq 



o \ sin(k - 1 )9 
q k - 1 



+ Constant 



(1+8) 



The equations of transformation corresponding to the 
other set of particular solutions from equations {\\) are 
obtained by replacing in equations (I4.8) the cosine by 
the sine and the sine by the negative cosine. 



The excluded cases k = 0 and k = ±1 are now 
treated. For k = 0, one set of particular solutions 
corresponds to a source and is 
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o o 



% = 9 



Equations (1|7) then yield 

P 



X = — - cos e 



y = ^ sm e. 

The other set of particular solutions corresponds to a 
vortex and is 

% = % 

Equations (Ij.7) then yield 

X = q- sm 9 
y = _ i cos e 

q 



For k = 1 with 

0^ = cos 9 
= - sin 9 

equations (l+y) yield 

^1 ^ r( -IT - ©1 cos 29 + ^ l( ^ — i + ^ I dq 
^ pq -1/' ^JV^ dq pq dq / 



3h 
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- 12. On ) sin 20 - U^i + ^1 6 

2 V q pq 



With the use of equations (Ip) and (5), 
1 - (1 - T)P*1 



X 



1 

1 ~ 2 



1 - 



(P + 1)t(1 - T)P 



COS 26 + log q 



2(3+1 



(1 - T)P 



- 1 



> (lf9) 



1 



1 - 



1 - (1 . T)^^-'- 

(p + 1)T(1 - r)^ 



sin 29-9 



where g(T) = g^(T) by equation (59) s^id is evaluated 
in part I (reference 1). 

For k = 1 with 

= P, sin 9 



\ = Q]_ cos 9 



_ 1 



1 - 



yi = 2 



1 - 



1 ~ (1 - T)P-^^ 
(p + 1)T{1 - t)' 

1 - (1 - T)P+1 



sin 29+6 



(p. + 1)T(1 - T) 



cos 29 + log q 



p+1^ ' 2P+1 



1> (50) 



(1 - T)P 



- 1 
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For k = -1 virith 



^-1 = P-l 



COS 6 



^Ki ^-1 sin e 



equations (ii.7) yield 

. ^ = 1 (Li . 
"•^ 4 \ q pq / 

4 \ q pq / 



cos 29 + ^- - 



po 

2jVq dq pq dq / 



1 -/I dP.i 



dq 



sin 28 - "1 



P. 



oq 



With the use of equations (l6) and (5)» 



2 



$p + 2 1 _ p 
P+l (l-T)P P+l 



(1 + PT) 



cos 2e 



log q + ^-^tS 



3(3-1-2 2p + 1 
8(P + 1) pa 2 



(1 - T)l 



- 1 



> (51) 



3P-I-2 1 p_ 

p+l {1-T)^^"P + 1 

1 



(1 + pt) 



sin 28 



ka 2 



e 
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For k = -1 with 



= - P.i sin 

= Q_.j^ COS 9 



2 



3P + 2 1 ^ 

p+1 (l-T)P (3 + 1 



(1+ PT) 



sin 29 



l+q^ Lp + 1 (1 -T)P p + 1 



(1 + PT) 



cos 29 



> (52) 



1 log q + ^ ^ g(T) 



+ 3P + 2 2p + 1 
8(P + 1) pa 2 



(1 - T) 



- 1 



Ringleb (reference 2) gives an example of the flov/ 
of a compressible fluid around a semi-infinite line. An 

exsraination of Ringleb »s stream function ^ = q sin 9 

shows that it is a linear combination of \|/-, and n ; 
that Is, ^ .'^ 



sin 9 
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In fact, all the external angle flows (l ^ -k < <» ) are 
nonuniquej for, in view of the discussion preceding 
equation (11), a general form of if_^ is 



*-k ^ 



At Y^(t) + Y_^Jt) I sin ke 



where A is en arbitrary constant. 



Observations on Limit Lines 

In the present section there are reviewed briefly 
certain conditions, discussed by Tollmlen (reference 9) 
and Rlngleb (reference 10), with regard to possible limi- 
tations on the potential flow of an adiabatic compres- 
sible fluid. 

. Consider the family of streamlines 



\l/(6jq) - Constant 



Then along a streamline 



d\|; = |i de ^ 1^ dq = 0 

C9 OQ 



and, from equations ikl):^ the line elements along a 
streamline ere 



q 



cos 6 
6^^l/^B 



dq 



dq 



(53) 



Singular points along a streamline are characterized by 
the vanishing of the common factor of equations (53): 
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(0 



d\!/ 



and 



d\l/ 



vanish, . the 



and the source for which ^ - 



(Stagnation points at which 
vortex for which. 5^ " ^ 

are excluded from this discussion. ) Observe now from 
equations (1|7) that the Jacobian of the transformation 
from the hodograph variables 8 and q to the physical 
plane variables x and y is given by 



fliT\ = Po^ fd£ d± d£ 6±\ 
V0.q>^ pq2 v59 dq " 5q 59/ 



q V 0 / 



0 



(55) 



Thus, the vanishing of the Jacobian is equivalent to the 
condition for the existence of a singular locus for the 
family of streamlines 



^(6,q) = Constant 

This singular locus consists of points at which the 
streamlines undergo an abrupt change of curvature and 

means, physically, that the acceleration q 4^ of a 
fluid particle is infinite at such points* 

Both Ringleb and Tollmien. have shown that the 
singular locus for the streamlines is also the envelope 
of the Mach lines in the plane of flow. The Mach lines 
are related to the streamlines in such a way that the 
component of the fluid velocity normal to a Mach line is 
equal to thie local velocity of sound. The Mach lines 
are identical with the so-called characteristic curves 
of the second-order partial differential equations 
for 0 and xl/ and are the integral curves of the ordi- 
nary differential equation 
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d8 



- \{lt - l) dq2 = 



or 



d9 = ±|-v^^ - 1 dq 



(56) 



The real solutions o.f this differential- equation inter- 
preted in the physical x, y plane yield the Mach lines 
for a given flow. The solution of equation (56) is 



9 - = ± 



(57) 



where i"' - 



+ 1 



and vvhere 9^ ass'oraes the values 



of 9 along ths M = 1 line for a gp.ven flow. 
It is recalled that the function 

H = log q + h(T) 



introduced in part I (reference 1) in connection with 
the geometric -mean tj-po of velocity correction formula, 
is a solutlori of the dlf f erer.tlal equation 



dH 



/ ? 



dq 



in the subaonic renge. Observe that a continuation of 
the function H into the supersonic range is given by 
equation (36) as 



d9 = ± ^ 



dq 



In the supersonic range, the function H = 9 - 9^ can 
thus be intei^jreted as the hodograph of the Mach lines 
for a given flow. 



1^0 
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The differential line elements dx and dy for the 
Mach lines in the physical plane are now given. Prom 
equations (J+7) and (56), the line elements along a Mach 
line for a given flow are 



dx = cos 9-sin e)(|^lU'2.i ||^a, 

dy « ^(t^^pTl sin e +COS e^(^tly^~l H^dq 



(58) 



3ingxj.lar points along a Mach line are characterized by 
the vanishing of the common factor of equations (58) 



dq 



(59) 



Equation (59 ) represents in the plane of flow two pos- 
sible singular loci or ^^lirnit lines'' for the tv/o families 
of Mach lines associated with the plus end minus signs 
in equation (56)* Clearly, the two singular loci cannot 
occur simultaneously since the two conditions cannot be 
satisfied simultaneously. Observe that equation (59) is 
equivalent to the vanishing of the Jacobian given in 
equations (55)- Thus^ the vanishing of the Jacobian is 
not only the condition for the existence of a singular 
(cusp) locus for the streamlines but also the condition 
for the existence of a limit line (envelope) for the 
Mach lines. 

The existence of a singular locus may be locked upon 
as being equivalent to the vanishing along a curve of 

the Jacobian the transformation from the 

hodograph-plene variables 9 and q to the physical- 
plane variables x and y. It is remarked that 

singular solutions exist for which the Jacobian j(^^^ 

Vx,y/ 

of the transformation from the physical -plane variables x 
and y to the hodograph-plane variables 6 and q 
vanishes identically in a region of the physical plane. 
In this case, as Tollmien pointed out, 9 and q are 



NAG A ARR No, ri;l29 



no longer independent variables and the flow cannot be 
described In the hodograph plane* Examples of these 
"raissed flows" are the solutions of Meyer (reference 11) 
for supersonic flow inside and outside sharp angles. 

It is of special interest to apply the condition 
for the vanishing of the Jacob i an to the particular 
solutions 0^ and treated in the early part of 

this paper. Th.e expression for the Jacobian for a 
particular solution 

0. = cos ke 



is, with the use of equations (5), for k 7^ 0, 



2 Vd8 dq 



6^ dii\ 
dq 69/ 



Pk^ ^ 



2 

in^e + (l-M^) Q^^ CO 



s^ke 



(60) 



Clearly, this expression for J is positive in the sub- 
sonic range M < 1. At the sonic value M = 1, ^ 0 
(see table following equation (25)) and J is again 
positive. At the first zero of P^^ in the supersonic 
range M > 1, \ ^ ^> hence, J is negative. The 
values of M, for all the pairs of values 6, M for 
which the Jacobian J vanishes, therefore lie betv/een 
M = 1 and the value of M at the first zero of F^^ 

(or S^^) in the supersonic range. 

By means of the relation 



Pk = f w 



i+2 
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the vanishing of the Jacobian yields 



cot k6 = + 



Sk 



(61) 




Squation (6l) is the relation for pairs of values 9, M, 
which interpreted in the physical x, y plane constitute 
the limit line for the particular flow {2^, \|/^, The 

values of M that satisfy equation (6l) accordingly lie 
between M = 1 and the value of M at the first zero 
of in the supersonic range. 

This paper is closed with the following remarks on 
limiting values of M in connection with the use of 
velocity correction formulas. The limiting local values 
of M in the case of uniform flow past a prescribed 
boundary, in general, depend on shape parameters. The 
use of a velocity correctlo'^ formula, however, yields a 
constant limiting valu6 of M that depends only on the 
particular correction formula Ui:$fe.d.? The geometric-mean 
correction formula yields the- value M = 1; the approxi- 
mation of Temple and Yarwbod yields M = 1.55 J ^nd the 
arithmetic -me an correction formula given in part I 
(reference 1), which is based on a linear combination of 
a source (limiting value M = 1) and a vortex (limiting 
value M = oo) or a spiral flow, yields the value M = 1.3 
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